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LATITUDINAL VARIATIONS IN THE 
NEUTRAL ATMOSPHERIC DENSITY 


George P. Newton and David T. Pelz 


ABSTRACT 


In-situ atmospheric density measurements from the Explorer 
32 density gauges for the time period May through October 1966 
reveal latitudinal density structure in the neutral thermosphere 
during geomagnetically undisturbed times. At altitudes 'between 
400 and 700 km the densities at northern geographic latitudes 
(>55°) are at least 1.5 times the equatorial densities during the 
daytime, and greater at night. The northern auroral zone atmos- 
pheric densities have been observed to be as much as a factor of 
five greater than the equatorial densities at altitudes between 450 
and 600 kilometers, even on geomagnetically quiet days. These 
results suggest that auroral zone and high latitude heating is a 
permanent feature of the summer thermosphere. 
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LATITUDINAL VARIATIONS IN THE 
NEUTRAL ATMOSPHERIC DENSITY 


INTRODUCTION 

Early conclusions regarding the latitudinal variation of the neutral thermos- 
pheric density were based on atmospheric densities derived from variations in 
the orbital periods of satellites. These conclusions and conclusions from pre- 
vious density gauge measurement results are presented and summarized in the 
following paragraphs . 

Lidov (1958) reported that near 225 km altitude the northern values of 

a/h (product of density and square root of the density scale height)- were roughly 
1.5 times the southern values. Schilling and Whitney (1959) reported data ob- 
tained during July to November 1958 in the altitude range 257 to 270 km, and 
concluded that there was no appreciable variation of average density between 
latitudes 49° N and S. At about 50° N and S latitudes however, they observed a 
sharp discontinuity in the density with the northern hemisphere value about twice 
the southern hemisphere value, and assuming the phenomenon to be real inter- 
preted the result to be a variation of density with latitude or time. Champion 
and Minzner (1959) reported preliminary results which suggested that the density 
in the altitude range 170 to 230 km may be lower in the summer than in the 
winter. Groves (1961), for altitudes below 300 km, concluded that the atmos- 
pheric density varied less than 20% with latitude except for a possible 60% in- 
crease in density near the winter pole. King-Hele and Walker (1961) concluded 
that there was no systematic latitudinal density variation exceeding a factor of 
1.5 at heights between 200 and 600 km between latitudes 50° N to 50° S. Paetzold 
and Zschomer (1961) stated that a weak dependency of density on latitude and 
season exists such that in summer at 200 km height the polar region density is 
15% higher than at lower latitudes , while in winter the density decreases by 10% 
towards the pole. May (1963) reported that the diurnal density variation at 205 
km altitude was smaller at higher latitudes. In a later paper (May, 1964) he - 
concluded that there was a decrease in density of approximately 30% with in- 
creasing latitudes at 1400 hours local time for altitudes below 300 km. Roemer 
(1966) found a latitude-dependent seasonal density variation of ±25% at 39° lati- 
tude and 690 km altitude with the greater density in the winter hemisphere. More 
recent findings (Keating and Prior, 1967), (Jacchia and Slowey, 1967) from satel- 
lite drag observations of satellites with perigees above 500 km, were that a win- 
ter helium bulge existed at these altitudes. Jacobs (1967) observed a daily peri- 
odicity in the accelerations of low perigee satellites (below 250 km altitude) 
which he attributed to an earth fixed density bulge at high latitudes. 

From Explorer 17 in-situ density gauge measurements, Newton et al. (1965) 
concluded that between 35° and 55° northern geographic latitude the average 
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latitude dependence of the atmospheric density has an amplitude of less than a 
factor of two. In a later paper (Newton, 1967) he concluded that the density at 
350 km altitude was a factor of two greater at 55 °N (summer 1963) than at the 
equator. 

If consideration of the differences in solar activity due to different times in 
the solar cycle corresponding to the observations are neglected , the earlier 
conclusions may be summarized as follows . There is general agreement that 
the latitudinal density variation for midlatitudes is less than a factor of two, al- 
though it is uncertain what the variation is. There is disagreement as to wheth- 
er the density at higher latitudes is greater than the equatorial density in sum- 
mer or in winter. The magnitudes of the density variations with latitude appear 
to be variable and may be a function of altitude. The reasons for the differences 
in the observations are as yet unexplained. 

Current atmospheric models describe the latitudinal density variation as 
that resulting essentially from a subsolar bulge with the density at a given alti- 
tude decreasing as the distance from the subsolar point increases. 

We present new direct measurements from the Explorer 32 density gauges 
which show increased density at high northern latitudes in May to October 1966 
contrary to the model atmosphere density behavior. This high latitude enhance- 
ment is believed to reflect an important source of heating of the high latitude 
and auroral zone thermosphere that may be a permanent feature at least in the 
summer hemisphere. 


DENSITY EXPERIMENT AND DATA ANALYSIS 

Direct, in-situ density measurements from the Explorer 32 atmospheric 
density experiment have been utilized in this analysis. Newton et al. (1968) 
have described the experiment and the density gauge response characteristics 
have also been reported (Pelz and Newton, 1967). 

The data we have considered are shown in Table 1 and were selected for 
days (1) which were geomagnetically quiet or undisturbed, and (2) for which den- 
sity measurements at different latitudes at common altitudes were obtained. 

The data cover the time period May through October 1966 and approximately one- 
half of the days meeting the above criteria were used in the analysis. However, 
the data reported are generally representative of the total data available during 
this time period. 

Figure 1 shows the measured atmospheric density versus altitude for 8 June 
1966, and is an example of a selected days data. It should be noted that turn-ons 
185 through 188 all occurred within 23 minutes of real-time, a time duration 
short compared to expected time lapses required for atmospheric density 
changes . 
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An individual display similar to Figure 1 was constructed for the analysis 
of the data corresponding to each day. For the altitude ranges in which both 
equatorial latitude data (between ±16° geographic latitudes) and higher northern 
latitude data existed, the ratio of northern latitude density to equatorial density 
at a selected altitude was calculated. This ratio was then multiplied by a density 
ratio computed from the Jacchia (1964) model (0° latitude, A p = 3, F10.7 = 100 
units) as follows: For the selected altitude the model density corresponding to 
the local time of the equatorial data was divided by the model density correspond- 
ing to the local time of the northern latitude data. This procedure adjusted the 
equatorial data to correspond to the local time of the northern latitude data. 


ACCURACY 

The accuracy of the density ratios is believed to be influenced primarily by 
the experiment precision. The agreement between the densities measured on 
turn-ons 188 and 193, as seen in Figure 1, is believed to be representative of 
the repeatability (>90%) of the density measurements. It should be noted that 
these turn-ons have a region of common altitudes, latitudes and local times. 
Neglecting any relative composition changes in the atmosphere as a function of 
latitude, the density ratio error is less than 20%. 

Density gauge measurement interpretation has some dependence on the at- 
mosphere composition. The density values in Table 1 and those shown in Figure 
1 were obtained by interpreting the data in terms of anatomic oxygen atmosphere, 
making appropriate recombination assumptions as previously discussed (Newton 
et al. 1965). The calculated possible impact of compositional corrections is 
summarized in Table 2. 

Appropriate composition corrections will be applied to the data at a later 
date when final analysis of the Explorer 32 and other mass spectrometer results 
has been completed. Errors in the density ratios due to incorrect atmospheric 
composition assumptions could be important for altitudes above 600 km, since 
this altitude may be near the transition level for atomic oxygen and helium and 
the transition altitude may vary with latitude. However, we believe the effect of 
composition corrections on the conclusions of this paper are small since most of 
the data considered are below 600 km. 


RESULTS 

The local time corrected density ratios have been 3 c pci r<xt cd into three ten 
degree wide latitude bands according to the geographic latitude of the northern 
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hemisphere data, and are shown in Figures 2-4 as a function of altitude. If the 
data are separated into comparable bands of geomagnetic latitude and plotted as 
a function of altitude, there is only a small change in the appearances of Figures 
2-4. Local times corresponding to the high latitude data are coded by symbols. 

It should be remembered that these data cover the time period May through 
October 1966, during which time the monthly average F 107 was 100 ±22 units, 
and correspond to geomagnetically undisturbed times. 

From Figure 2, it is seen that the low midlatitude density can be 20% great- 
er than the equatorial density during the day for altitudes between 300 km and 
400 km. 

The data in Figure 3 show that the density at midlatitudes is generally 
greater than the equatorial density over the altitude range 300 km to 500 km. 

The average density ratio is near 1.3 for these altitudes. 

The data in Figure 4 show the high northern latitude (55° N - 65° N geo- 
graphic) density to be greater than the equatorial density. Although during the 
day the average density ratio is near 1.4 to 1.5 and increases at night to be near 
a factor of two for altitudes between 370 km and 600 km, it can, as shown, be as 
large as a factor of 10. Generally, the higher density ratios (greater than two) 
correspond to locations near the auroral zone. In particular, all density ratios 
greater than three correspond to the high northern latitude measurements at 
locations in or near the auroral zone. In two instances these larger ratios cor- 
respond to the same geomagnetically quiet day. 


CONCLUSIONS AND DISCUSSIONS 

Our conclusions from the above analysis are: 

(1) There is a latitudinal variation greater than 30% in the neutral atmos- 
pheric density at altitudes between 300 km and 700 km for latitudes 
greater than 45° N geographic, with the summer northern latitude den- 
sities being greater than the equatorial densities. 

(2) The noi’thei’n high latitude (>55°) to equatorial latitude density ratio is 
smallest during the day (near 1.3) and increases at night to values 
greater than two. 

(3) Large density ratios (>3) are observed only near or in the auroral zone. 
These ratios are observed even on geomagnetically quiet days. 

These results are in agreement with those obtained from analysis of Ex- 
plorer 17 data and reported previously (Newton et al. , 1965, 1967). The early 
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preliminary Explorer 32 data analysis conclusions (Newton, 1967) were invali- 
dated by errors in the preliminary satellite world map which was used. Our 
conclusions are consistent with those of Schilling and Whitney (1959) and Paetzold 
and Zschomer (1961) in that the northern summer hemisphere density is greater 
than the equatorial density. Our observations are also in accord with the results 
of Jacobs (1967) , although the relationship between the densities at 500 km and 
200 km altitude at these latitudes is not clear. 

We have not yet determined the role, if any, of a seasonal variation in these 
results . We are currently analyzing additional Explorer 32 data covering the 
time period October 1966 to January 1967, in an attempt to determine the sea- 
sonal variation of the atmospheric density. It is clear, however, that the varia- 
tions observed and reported here are generally much larger in magnitude than 
and opposite in direction to the latitudinal behavior of the neutral thermosphere 
density predicted by current atmospheric models (Jacchia, 1965) (Friedman, 
1967). 

A possible explanation for the observed latitudinal density structure requires 
a general heating of the high latitude summer thermosphere even during geo- 
magnetically quiet times. This heating is supported by a preliminary analysis 
(to be reported later) of the density scale heights provided by the same experi- 
ment covering the same time period. Such heating is estimated to be a signif- 
icant fraction of the solar EUV heating on the basis of a comparison of the 
magnitude of the latitude density variation with the magnitude of the diurnal 
density variation. It may also be relevant to consider the latitudinal variation 
of the neutral thermosphere density as related to the observed light ion depletion 
and heavy ion abundance at high latitudes reported by Taylor, et al. , (1968). 
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Figure 1 . Atmospheric density versus altitude measured on June 8, 1966. Each line represents nearly one 
hundred density measurements. The turn-on number, latitude and local time ranges covered during the pass 
are shown beside the line. The start and end of each turn-on are represented by the circle and cross, 
respectively. The F iq j was 97 units and the Ap was 4 q. 
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Figure 2. Ratios of the 35°N-45°N geographic latitude densities to the 
equatorial densities versus the altitudes at which the ratios were determined. 
The ratios have been corrected for local time differences in the two density 
measurements. Local times corresponding to the higher latitude data are 
coded by symbols. 
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Table 1 


Table 1 is a listing of the Explorer 32 density measurements used to obtain 
the results presented in this paper. 

Col. 1 is the station code for the minitrack stations recording the data or 
the remote location identification in the case of tape recorded data. 

Col. 2 is the orbit number corresponding to the data. 

Col. 3 is the turn-on number corresponding to the data. 

Col. 4 is the date the data were obtained. 

Col. 5 is the average GMT corresponding to the data. 

Col. 6 is the daily F10.7 in units of 10' 22 W. /m 2 e/s corresponding to the date 
indicated. 

Col. 7 is the daily A p corresponding to the date indicated. 

Col. 8-12 are a listing of the values of the geographic latitude, longitude, 
local solar time, altitude and density respectively. 
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Table 2 gives the multiplicative correction factors to be applied to the densities 
given in Table 1 for several relative atmospheric compositions . 


Table 2 


Relative Atmospheric Composition 

Multiplicative Correction Factor for 
the Density Values in Table 1 
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